Minimal Single-Particle Hamiltonian for Charge Carriers in Epitaxial Graphene on 

4tf-SiC(0001) 
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We present a minimal but crucial microscopic theory for epitaxial graphene and graphene nanorib- 
bons on the 4//-SiC(0001) surface - protopypical materials to explore physical properties of graphene 
in large scale. Coarse-grained model Hamiltonians are constructed based on the atomic and elec- 
tronic structures of the systems from first-principles calculations. From the theory, we unambigu- 
ously uncover origins of several intriguing experimental observations such as broken-symmetry states 
around the Dirac points and new energy bands arising throughout the Brillouin zone, thereby es- 
tablishing the role of substates in modifying electronic properties of graphene. We also predict that 
armchair graphene nanoribbons on the surface have a single energy gap of 0.2 eV when their widths 
are over 15 nm, in sharp contrast to their usual family behavior. 

PACS numbers: 73.20.-r,81.05.Uw,68.35-p,71.20-b 



Graphene has attracted immense interests because of 
the unusual relativistic energy dispersions with the chiral 
massless Dirac fermions near the Fermi level (E F ) 
3]. The direct observation of such peculiar quasiparticle 
spectra in graphene is particularly important not only for 
understanding its novel physical properties 0, 0, 0, 0] but 
also for practical applications 0,11]. So, several measure- 
ments through the high resolution angle resolved pho- 
toemission spectroscopy (ARPES) are performed on a 
single layered epitaxial graphene lying on the silicon car- 
bide (0001) surfaces @0, 0, 0, [HI- Surprisingly, the 
reported quasiparticle spectra reveal anomalous energy 
dispersions around the Dirac energy point (Ed) indicat- 
ing highly renormalized bands 
Many experimental I8l.l9lll0l.lll 
and theoretical studies [13. |2 
address these important problems but there is still no 
consensus on origins of the anomalous spectrum. More- 
over, there is no clear understanding on other anomalous 
ARPES observations such as the broken six-fold sym- 
metry near the Ed and new distorted hexagonal energy 
bands around Brillouine zone (BZ) center 0,0,1,0. 

The complex interfacial structures arise when epitaxial 
graphene is grown by annealing SiC surfaces [271 ]. Dur- 
ing the thermal decomposition of the surfaces, a layer of 
carbon atoms, called the buffer layer, forms first with- 
out exhibiting the typical linear energy dispersion of 7T- 




states near the Ep 27 



, |29J, |30J. Then, on top of 
the buffer layer, the clean honeycomb lattice of carbon 
atoms grows [II [H [H, El HH . The lattice mismatch 
between the SiC(0001) surface, buffer layer and graphene 
gives rise to the large scale surface reconstruction with a 
periodicity of 6\/3 x 6\/3~R30 o (in short, 6i?3) with re- 
spect to the SiC(0001) surface unitcell, which is observed 
in the low energy electron diffraction (LEED) measure- 
ments 27 , 28l US EiH . The scanning tunnelling mi- 



croscopy (STM) image, however, indicates an approxi- 
mate 6x6 periodicity [ll|, Ei , IS 16[ of the reconstruc- 
tion. We note that such superperiodic rearrangements 
shall impose constraints on possible theoretical models 
to explain the anomalous electronic structures mentioned 
above. 

In this paper, we show that the interactions between 
epitaxial graphene and the reconstructed layer under- 
neath it are the main driving forces to several anomalous 
features observed in recent experiments @, 0, 0, 0, ES 
H El ES El ES ES E3, El on electronic properties of 
epitaxial graphene. From the simulated ARPES spectra 
on the systems, it is shown that the symmetry break- 
ing at the Ed and the new hexagonal bands throughout 
the two-dimensional BZ have the same origins. Based 
on the model, we predict that graphene nanoribbons or 
finite size fragments of graphene 0, EH, [28| show an 
homogeneous energy gap contrary to their family behav- 



32j. Moreover, by extending our microscopic model 



for monolayer epitaxial graphene to bilayer one, we iden- 
tify the effects of the buffer layer to its e nerg y spectrum 
and characteristic energy gap @, 0, IS ES E3] ■ Our com- 
putational results indicate that the interplays between 
geometries and electronic structures are pivotal in al- 
tering global energy bands of graphene, notwithstanding 



that the many-body interactions [2j, 25j, |26j are expected 



to play some roles in modifying the quasi-particle spec- 
trum near the Ed locally. 

From the ab-initio pseudopotential density functional 
method [22], [IS [13], we find that the carbon atoms 
in the buffer layer with superperiodic 6i?3 unitcell are 
split into lattice matched regions, where carbon atoms 
have cr-bonds to silicon atoms of the 4ff-SiC(0001) sur- 
face, and their boundaries without the cr-bonds [22| (Fig. 
1(a)). The carbon atoms at the boundaries of the lattice 
matched regions exhibit an approximate 6x6 domain 
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satisfying the geometric constriction imposed simultane- 
ously by both LEED and STM measurements. To ex- 
plore various initial conditions for buffer layer formation, 
we shift the initial atomic coordinates of the buffer layer 
(shown in Fig. 1(b)) on top of the 4ff-SiC(0001) sur- 
face by either iai (Fig. 1(c)) or ia 2 (Fig. 1(d)) where 
a.i(2) is an unit vector of graphene. We find that the fi- 
nal relaxed atomic geometries for all initial coordinates 
are essentially same to each other except for minor differ- 
ences in the connectivity of 7r-electrons along quasi-6 x 6 
domain boundaries (Figs. l(b)-(d)). 

We build up a minimal (coarse-grained) microscopic 
model for interactions between epitaxial graphene and 
the buffer layer based on atomic and electronic struc- 
tures obtained from our first-principles calculation (See 
detailed method in [lij]). From the first-principles calcu- 
lations, it is found that the interactions between 7r-orbital 
states at the domain boundary of the buffer and ones 
in graphene play the most significant role to determine 
the electronic structures of the system while the states 
of atoms inside the domain and those under the buffer 
layer have negligible contributions to the electronic struc- 
tures near the Ef and Ed- Hence, it is sufficient to ap- 
proximate the buffer layer to the coarse-grained atomic 
configuration of the quasi-6 x 6 periodic connections of 
7r-electrons only (thick black lines in Figs. l(b)-(d)). 
Our tight-binding Hamiltonian for monolayer epitaxial 
graphene on coarse-grained buffer layer model (Fig. lb) 
can be written as 

n = -t G 53 4 c 3 ~ v gY1 4 c * - *b 53 h \ h ™ 

% (l,m) 
-V B J2 b lh-jJ2 C tbm + (c.C.), (1) 

where t G (2.70 eV) and ie (1-50 eV) are the near- 
est neighbour hopping amplitude between carbon atoms 
in graphene and those in the buffer layer respectively. 
The Ci and bi are annihilation operators for electron 
in graphene and buffer layer respectively. 7 (0.30 eV) 
denotes the interlayer hopping amplitude between the 
nearest neighbour carbon atoms belong to graphene and 
the buffer layer with the Bernal type stacking respec- 
tively. V G (0.35 eV) and V B (0.34 eV) describe the 
potential for graphene and the buffer layer considering 
charge redistributions due to the polar SiC surface. Since 
our model Hamiltonians are described within the single- 
orbital tight-binding approximation, we can readily ex- 
tend our model to explore the electronic structures of 
finite-sized monolayer and multilayer epitaxial graphene 
respectively which are beyond the reach of the first- 
principles calculations. 

We have found that the nearest neighbour inter- and 
intra-layer interactions between 7r-electrons in the coarse- 
grained atomic model (Eq. (1)) are sufficient to repro- 
duce the electronic energy bands obtained from the first- 



phene (b) *. 




3 z 




(d) M X^A 




FIG. 1: (color online) (a) Ball and stick model for fully relaxed 
atomic configuration of epitaxial graphene, buffer layer and 
4i/-SiC(0001) surface with 6R3 supercell. Epitaxial graphene 
(black) is located on top of the buffer carbon atoms (black and 
grey). The buffer carbon atoms consisting 6x6 domain are 
denoted by grey color and its boundary by black. The silicon 
and carbon atoms in the SiC are denoted by red and green re- 
spectively. (b)-(d) Coarse-grained minimal microscopic mod- 
els. The red thick arrows indicate the 6x6 supercell unit 
vectors and the blue ones 6-R3 unit vectors. Black lines repre- 
sent the connectivity of 7r-electrons of the buffer carbon atoms 
with approximate 6x6 domain boundaries and grey lines de- 
note the hexagonal network of epitaxial graphene on top of 
the buffer. 



principles calculations (Fig. 2) (for detailed comparisons, 
see [34|). We notice that inclusion of the nearest neigh- 
bour interlayer interaction between graphene and the 
coarse-grained buffer layer model already breaks the sym- 
metry between two sublattices of graphene. The simu- 
lated ARPES spectrum for monolayer epitaxial graphene 
agrees well with both our previous results _[22j from the 
first-principles calculation and ones @, 0, Hl^j E3] from 
experiments (Fig . 2) (for detailed ARPES simulation 
method, see [34]]). When the initial conditions for the 
buffer layer formation are varied and the resulting ge- 
ometries for the approximated 6x6 domain are slightly 
altered as shown in Figs. l(b)-(d), the simulated ARPES 
spectrum display essentially same structures with energy 
gaps at the Ed's and the midgap states as shown in Figs 
2(a)-(c). Hence, We confirm from the simple microscopic 
model that the sublattice symmetry-breaking interlayer 
interaction indeed opens a gap of ~0.20 eV at the Ed and 
the presence of midgap states give rise to high ARPES 
intensities inside the energy gap as well as the level re- 
pulsion between upper and lower Dirac cones (Fig. 2(d)). 

From simulated ARPES spectrum for bilayer epitaxial 
graphene, we find that the Ed approaches to the Ef 
and the energy gap decreases from 0.20 to 0.12 eV (Fi, 
2(e)), agreeing well with experimental observations [7|,ll 
17J Due to the charge transfer between the SiC surfaces 
and graphene, perpendicular electric fields exist on the 
multilayer epitaxial graphene opening an energy gap 35, 
36] . Though the fingerprints of midgap states are hardly 
visible for the bilayer epitaxial graphene, it is shown that 
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FIG. 2: (color online) Simulated ARPES spectrum near the 
Ed for various atomic models are drawn along the arrow 
shown in the inset of (a). The superimposed dotted lines 
are an ideal energy spectrum of graphene. Each simulated 
ARPES spectra in (a)-(c) corresponds to the geometry shown 
in Fig. 1(b)- (d), respectively. The positions of Ed and size 
of energy gaps are (a) -0.35 eV, 0.20 eV, (b) -0.33 eV, 0.19 
eV, and (c) —0.33 eV, 0.23 eV , respectively, (d) Schematic 
energy dispersion for epitaxial graphene with gap and midgap 
state. The straight lines are linear dispersion relations of ideal 
graphene. (e) Simulated ARPES spectra of bilayer epitaxial 
graphene near the Ed- Superimposed red and blue lines rep- 
resent energy bands of bilayer graphene with and without 
buffer layer respectively. The position of Ed and energy gap 
are given by -0.26 eV and 0.12 eV respectively. 



(a) 


J» '"*L 


1 


r 1 




^ & 






E F -0.8eV 


(c) 


< 




> 






E F 



(b) 












E D (E F -0.35 eV) 


(d) 


J w 






E F -0.35 eV 



FIG. 3: (color online) Simulated ARPES intensity maps for 
constant energy taken on monolayer epitaxial graphene at (a) 
-0.8 eV (b) -0.35 eV (E D ) and (c) 0.0 eV (E F ) respectively. 
The first Brillouine zone of graphene is denoted by dotted 
line in (a). The insets in (a) and (b) show the stereographic 
plots near the Ed ■ The small arrows in (b) indicate the second 
shortest reciprocal lattice vectors of 6i?3 supercell, which con- 
nect the AT-point to six faint replicas, (d) Simulated ARPES 
intensity maps taken on buffer layer at —0.35 eV (Ed) with- 
out graphene. The bar on the right side denotes the color scale 
of relative intensities from zero (white) to maximum (red). 



the interaction between the buffer and graphene flattens 
the characteristic Mexican-hat-shaped band j3f| at the 
bottom of the upper Dirac cone [Fig. 2(e) and Fig. IS 
(c) in 34]] and shifts the second subband in the lower 
cone downward by 20 meV (Fig. 2(e)). 

By simulating the two-dimensional ARPES intensity 
maps with various fixed energies (Fig. 3), we find that 
the interactions between the buffer layer and graphene 
indeed induce the symmetry breaking near the Ed and 
produce new hexagonal energy bands throughout the first 
BZ of graphene. Our simulated fixed-energy intensity 
patterns give a good agreement with experimental re- 
sults [1, UH, 0] (Fig. 3). Due to the two-source interfer- 
ence between photo-excited electrons from two equivalent 
atomic sites (d la Young's double slit) or helical nature of 
the carriers in graphene, the intensity patterns at the K- 
points show the typical crescent shape anisotropy [Tij ]. 
With the characteristic high intensity at the iiT-points, 
there exist six equivalent faint replicas around each K- 
point forming a smaller hexagon (inset in Fig. 3(a)). 
When approaching Ed, the highest intensities at the 
if-points become isotropic and three of six replicas be- 
come weaker compared to other three points (inset in 
Fig. 3(b)). The intensities of stronger three satellites 



amount to 7% of the main peak at the _ftT-points while 
those of the weaker to 2% agreeing with experiment re- 
sults qualitatively 0,H,0]. The replicas are connected by 
the second-shortest reciprocal lattice vectors of the 6i?3 
supercell (Fig. 3(b)) and the area of the smaller hexagon 
nearby is ^ x ^ of the first BZ of graphene 0] • When 
the interlayer interaction is set to zero intentionally in 
our simulation, we cannot find any symmetry breaking 
phenomena. Thus, we conclude that apparent six-fold 
symmetry breaking near the -ftT-points is due to the inter- 
action between the buffer and graphene. Together with 
six faint replicas, there are global faint features through- 
out the first BZ (Fig. 3(a)-(c)). We also find that the 
larger hexagonal structure around the T-point observed 
in the experiments 

UMM (25% of total area) originates 
from the buffer layer. By simulating constant energy map 
(Fig. 3(d)) and ARPES spectrum (Fig. 2S in 34]) of our 
coarse-grained buffer layer model, we show that the un- 
derlying faint features are reminiscent of the interlayer 
interactions between the buffer and graphene. 

Usually, the surface of epitaxial grap hene exhibits the 
finite-sized terrace patterns [H, [lil l28l|. Thus the quan- 
tum confinement effect [32] may play a role in deter- 
mining the energy gap [9(. To illustrate such an effect, 
we calculate electronic structures of epitaxial graphene 
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FIG. 4: (color online) (a) Atomic model for EAGNR of width 
W (here, 24-EANGR is drawn) on top of the buffer layer. The 
lines follow the convention introduced in Fig. 1. (b) Simu- 
lated ARPES spectrum for nanoribbons with W = 19.9 nm 
(163-EAGNR). The superimposed lines represent the energy 
band of 163-EAGNR having an energy gap of 0.72 eV without 
the buffer, (c) The energy gaps of 7V-EAGNRs as a function 
of width (N = 3p, 3p + 1, 3p + 2,p — positive integer). The 
dotted lines correspond to energy energy gaps belong to each 
family without buffer layer. 



nanoribbons with armchair edges on both sides (in short, 
EAGNR) on top of the buffer layer. Following the con- 
vention [32|, three families of EAGNRs are denoted by 
the number of dimer lines, N, i.e., N-EAGNR (N = 3p, 
3p + 1, 3p + 2 families, p is a positive integer). By using 
the same model for the single layer epitaxial graphene 
(Fig. 4(a)), we show that the EAGNRs have a similar 
ARPES spectrum to epitaxial graphene for energy gaps, 
midgap states, and level repulsion between upper and 
lower cones, respectively (Fig. 4(b)). Moreover, when 
the width of EAGNR is over 15 nm, the typical family 
behaviour of energy gaps [32| disappears completely and 
converges to the energy gap (0.20 eV) of two-dimensional 
epitaxial graphene (Fig. 4(c)). Hence, we can con- 
clude that the terrace patterns or finite-sized epitaxial 
graphene exhibit essentially the same electronic struc- 
tures of the ideal two dimensional one. On the other 
hand, the present calculation results indicate that epi- 
taxially grown graphene nanoribbons on the SiC surface 
will have the homogenous energy gap if the width is over 
15 nm. 

In summary, we have constructed the microscopic the- 
ory for expitaxial graphene on 47J-SiC(0001) surface in- 
corporating interactions between graphene and the sur- 
face. The simulated experimental observations based on 
the theory have been shown to explain the several atypi- 
cal aspects of epitaxial graphene from a single and unified 
view and, thus, shed light on understanding the quasipar- 
ticle spectrum of graphene in various circumstances. 

S. K and J. I acknowledge the support of the KOSEF 
through the SRC program (Center for Nanotubes and 
Nanostructured Composites). H. J. C. was supported by 
the KRF (KRF-2007-314-C00075) and by the KOSEF 
Grant No. R01-2007-000-20922-0. Y.-W. S. was sup- 
ported by Quantum Metamaterials Research Center No. 
Rl 1-2008-053-01002-0 and Nano R&D program 2008- 
03670 through the KOSEF funded by the Korean gov- 
ernment (MEST). Computational resources have been 



provided by KISTI (KSC-2008-S02-0004) and the KIAS 
Linux Cluster System. 



SUPPLEMENTARY INFORMATION 

Construction of minimal model Hamiltonian within 
tight-binding approximations 

In first-principles calculations, we expand the wave 
function with localized basis sets [H, [33| to handle a 
large number of atoms in the system (typically > 1600 
atoms). A single-^ for hydrogen, a singe-C plus polariza- 
tion for silicon and a mixed basis set with a single and 
double-^ for s- and p-orbitals of carbon atom have been 
used, respectively [22J. The Kleinman-Bylander's fully 
separable nonlocal projectors |37l| are used in the norm- 
conserving pseudopotentials [38] and the local density ap- 
proximation [39l | is employed in setting up the exchange- 
correlation potential. We have thoroughly tested our ba- 
sis set with other calculation parameters to reproduce 
the atomic and electronic structures of SiC, graphene, 
other model for epitaxial graphene studied in previous 
literatures [2^, [2l|, [23|, respectively. We modeled the 
4_ff-SiC(0001) substrate in the simulation with four al- 
ternating silicon and carbon atomic layers. Hydrogen 
atoms are introduced to passivate the dangling bonds in 
bottom of the slab. On top of the Si-terminated surface 
of 4£L-SiC(0001), one, two, and three graphene layers are 
placed for the buffer layer, monolayer graphene, and bi- 
layer graphene, respectively. 

Based on the low energy electron diffraction (LEED) 
experiments 0, [2^, [3l|, the large supercell with the 
6-\/3 x 6\/3i?30 o (in short 6i?3) periodicity (equivalent 
to 13 x 13 times graphene unit cell) is imposed to the 
calculations. The atomic positions are determined by to- 
tal energy minimization calculations until the forces on 
each atom are less than 0.06 eV/A while atoms belonging 
to the last two silicon and carbon layers are fixed to the 
bulk atomic structure of 4iJ-SiC. In the electronic struc- 
ture calculations after geometric optimization process, we 
use 2x2 k-point sampling in lateral directions and set 
very large size of vacuum (50 A) in surface normal direc- 
tion to prevent the spurious dipolar interactions between 
terminated slab geometries in supercell configuration. 

During the high temperature thermal decomposition 
process 



12, US I2a 13J, the detailed atomic 
structures and shape of the domain may vary from sam- 
ple to sample. However, the domain boundary of 7r- 
orbitals of the buffer carbon atoms and the exact 6i?3 
periodicity remains the same as observed in many exper- 
iments 0> E3, 3 IS 31 1 and shown in our calcu- 
lations [22 1 . Hence, it is sufficient to approximate the 
buffer layer to the coarse-grained atomic configuration of 
the quasi-6 x 6 periodic connections of 7r-electrons only. 
We also found from the first-principles calculations that 
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the interactions between 7r-orbital states at the domain 
boundary of the buffer and ones in graphene play the 
most significant role to determine the electronic struc- 
tures of the system (22j. It is also noticeable that elec- 
tronic states of atoms inside the domain and those under 
the buffer layer have negligible contributions. 

Hence, based on results from the first-principles cal- 
culations, our tight-binding Hamiltonian for monolayer 
epitaxial graphene on coarse-grained buffer layer model 
(Fig. 1(b)- (d)) can be written as 

i (l,m) 

-Vb$>|&,-7E 4h m + (c.c), (2) 

where £g(= 2.70eV) and £b(= 1.50eV) are the near- 
est neighbour hopping amplitude between carbon atoms 
in graphene and those in the buffer layer respectively. 
The Ci and bi are annihilation operators for electron in 
graphene and buffer layer respectively. j(= 0.30eV) de- 
notes the interlayer hopping amplitude between the near- 
est neighbour carbon atoms belong to graphene and the 
buffer layer with the Bernal type stacking respectively. 
V G (= 0.35eV) and Vb(= 0.34eV) describe the potential 
for graphene and the buffer layer considering charge re- 
distributions due to the polar SiC surface. These pa- 
rameters were found to be enough for fitting our first- 
principles energy bands of the buffer layer and monolayer 
epitaxial graphene, respectively (Fig. IS (a) and (b)) 



while the second and third nearest neighbour interlayer 
interaction terms improve the agreements a little. 

For bilayer epitaxial graphene, we introduce another 
interlayer hopping amplitude A between the nearest 
neighbours carbon atoms belong to each graphene layer 
(Bernal type stacking) and potential shifts for each 
graphene layer, Vq and Vq with respect to the Fermi 
level. The Hamiltonian for the interaction is written as 
H' = Tto + Hi, where 

Hi =-t G J2 d ldj-V G J2 d Ui- A E c l d J + ( c - c -)- (3) 

Here, di is the annihilation operator for electron in second 
graphene. We fit the energy spectrum of bilayer epitaxial 
graphene obtained by our model Hamiltonian to the first- 
principles calculation results (Fig. IS (c)) and found A = 
0.35eV being quit e similar to one of g raphite, but smaller 
than 0.48 eV [17| and 0.46 eV [18] obtained in recent 
experiments. 

As shown in Fig. IS, the agreements between energy 
spectrums obtained by the first-principles calculations 
and ones by our model Hamiltonians are excellent. The 
dense flat bands near the Fermi level in Fig. IS from the 
first-principles calculations are found to originate from 
localized states inside approximate 6x6 domains of the 
buffer layer and localized states in 47J-SiC(0001) sur- 
faces underneath the buffer. We found that such local- 
ized states forming flat bands are buried deep inside the 
surfaces and do not contribute to the simulated ARPES 
spectrum. 




FIG. 5: Comparison between electronic energy bands from 
the first-principles calculation and tight-binding approxima- 
tion. Electronic energy bands of (a) the buffer layer on top of 
4_ff-SiC(0001) slab model, (b) monolayer epitaxial graphene 
and (c) bilayer epitaxial graphene on top of the buffer layer 
with 4/f-SiC(0001) slab, respectively. The spectrum are 
drawn along the high symmetric lines of the first Brillouine 
zone of 6-R3 supercell. The dotted red lines are obtained by 
the first-principles calculations on the atomic models includ- 
ing one (two) graphene layer, buffer layer, and 4ff-SiC(0001) 
slab. The solid black lines are obtained by our minimal micro- 
scopic model Hamiltonian within the tight-binding approxi- 
mations (Eq. 1). It is found that the dense flat bands near 
the Fermi levels are from the states isolated inside the ap- 
proximate 6x6 domains of the buffer layer and SiC surface 
underneath it, which do not contribute to ARPES intensities. 



Simulation of angle resolved photoemission 
spectroscopy (ARPES) intensities 

We use the Fermi golden rule to simulate ARPES 
spectra 0, |4(J. The transition probability (I) from 
an initial Bloch state (^i) to an outgoing electron state 
(\f r / = e lp r ) is written as 

I ~ K*/|Wph_«i|*i}| 2 a(fiw + Ei - E f ) (4) 

where 7i p h-ci is the photon-electron interaction Hamilto- 
nian, uj is the frequency of incident photon, and £";(/) is 
an energy of ^i(f)- Tight-binding wavefunction obtained 
for the system is given by, 
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a,-(k) 



T)e 



ik-T 



(5) 



where cy(k) is an amplitude for the 7r-orbital located at 
Xj in the 6\/3 x 6\/3R30 o supercell with unit vector T 
and k is the crystal momentum of electron. 

By using the dipole approximation, 7i p h-ci — Ae lq ' x , 
the transition probability will be expressed as 

2 



I(P) 



cij(k)e" 



5(Huj + Ei- E f ) (6) 
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/f(A" 1 ) k(A- 1 ) k(A-i) A- (A- 1 ) 

FIG. 6: Simulated energy spectrum of epitaxial graphene and buffer layer along various directions in the first Brillouine zone 
of graphene. (a) Simulated ARPES intensity map taken on monolayer epitaxial graphene at the Dirac energy point {Ed)- 
(b) From left to right panels, simulated ARPES spectrum along the g\, g^ and g$ shown in (a) respectively. Together with 
strong ARPES intensities from graphene, there are several dispersions originating from the buffer layer, (c) Simulated ARPES 
intensity map taken on minimal buffer layer model at the Ed- (d) From left to right panels, simulated ARPES spectrum along 
the 61, &2 and 63 shown in (c) respectively. It is noticeable that all characteristic ARPES spectrum of the buffer shown in each 
panel of (d) appears faintly in the corresponding one of (b), respectively. 



where p is the momentum of outgoing electron and cj) p is [5] 



the Fourier transformation of atomic orbital 4>(r), defined 
by 4>p — J 4>ji r ) e ~ lp r dr (set by constant). The crystal 
momentum k belong to the first Brillouin zone of the 
6-\/3 x 6\/3i?30 o supercell and satisfys momentum con- 
servation condition in the surface parallel direction (i.e., 
k+ G = p|| for the reciprocal vector G) of the supercell. 
We include the attenuation factor for the ARPES inten- 
sity considering photoelectron mean free path (~5A) to 
the surface normal direction. The conservation of en- 
ergy imposed by the 5-function is replaced by Lorenzian 
function, ~ (hui+E -E ■)' i +v' i > with a broadening (T) of 30 
meV. We simulate the ARPES spectrum (shown in Figs. 
2 and 4(b)) with a photon energy of 50 eV and the two- 
dimensional constant energy maps (shown in Fig. 3 and 
Fig. 2S) with a photon energy of 100 eV @, 0). 
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